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1.  IHTRODUCTION 

The  x^enomenon  of  cavitation  was  described  by  Euler  (l)*  in  175^  in 
his  theory  of  turbines.  He  noted  that  an  insufficient  pressure,  or  even 
a  large  pressure  in  a  perfect  fluid,  can  cause  disagreement  between  theory 
and  experiment  and  can  result  in  a  state  of  zero  resistance.  Since  Euler's 
introduction  to  the  concept  of  the  formation  of  cavities  in  a  fluid,  the 
phenomenon  has  been  the  subject  of  a  largs  rumber  of  scientific  investl- 
gatlons(2) . 

Cavitation  is  a  dynamic  process  involving  the  formation  and  collapse 
of  vapor-filled  cavities  in  a  liquid  in  motion.  In  normal  liquids  these 
cavities  form  if  the  local  pressure  drops  below  the  vapor  pressure.  Con¬ 
versely,  a  collapse  occurs  when  a  cavity  is  transported  into  a  region  of 
hi^er  pressure.  The  collapse  process  can  be  strongly  influenced  by  the 
collision  between  a  cavity  and  a  rigid  object. 

Essentially, there  are  two  types  of  cavitation.  In  the  first  the 
cavity  remains  fixed  with  respect  to  a  bounding  surface  and  usually  is 
observed  in  hydraulic  machinery.  As  Euler(l)  and,  later  BaUTiaby(5)  noted,  an 
important  consequence  of  this  is  a  decline  in  machinery  efficiency.  In  fixed 
cavitation  the  main  flow  leaves  the  guiding  surface  and  follows  a  free  tra¬ 
jectory  which  usually  returns  to  the  surface  at  some  downstream  point.  The 
fixed  cavity  occupies  the  space  between  the  solid  guiding  surface  and  the  free 
liquid  surface.  Knapp(4)  has  suggested  that  this  type  of  cavitation  is  quite 
different  from  the  second  type,  which  is  characterized  by  discreet  cavities 


*  Numbers  in  parentheses  refer  to  references  at  the  end  of  the  paper 
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moving  with  the  fluid.  The  model  for  this  type  of  cavitation  is  the 
collapse  of  a  spherical  huhble. 

Even  though  a  typical  cavity  in  water  is  approximately  only  one 
centimeter  in  diameter^  it  may  collapse  violently  and  generate  pressures 
of  the  order  of  hundreds  of  kilohars.  It  is  not  clear  which  aspect  of  the 
collapse  is  responsible  for  damage^  such  as  the  pitting  of  propellers. 

After  the  collapse  a  pressure  wave  is  propagated  outward  in  the  form 
of  a  shock  wave,  ^Ich  may  cause  some  damage  to  neighboring  rigid  sur- 
faces.  On  the  other  hand,  the  work  of  Komfeld  etnd  Suvorov(5)  indicates 
that  the  real  damage  is  done  when  the  cavity  strikes  the  object  and  not 

by  a  pressure  wave  initiated  by  a  collapse. 

« 

To  investigate  the  phenomenoncf  a  cavity  collapse  completely,  account 
should  be  tedcen  of  the  vapor  inside  the  cavity  and  the  effects  of  viscosity, 
surface  tension,  and  compressibility.  The  stability  of  the  process  should 
also  be  considered.  So  far,  theoretical  investigations  have  treated  these 
properties  in  isolation  (or  at  best  have  considered  two  simultaneously, 
neither  of  which  was  compressibility)  and  the  present  work  is  an  extension 
of  considerations  which  have  only  dealt  with  liquid  conpresslblllty. 

One  of  the  earlier  theoretical  investigations  concerning  the  collapse 
of  a  spherical  cavity  was  made  in  1917  by  Lord  Rayleigh(6).  He  refers  to 
and  extends  the  work  of  Besant(7}  who  calculated  the  pressure  distribution 
at  the  time  of  collapse.  Eesant's  statement  of  the  problem  is  .  •  .  "An 
infinite  mass  of  homogeneous  incompressible  fluid  is  acted  upon  by  no 
forces  at  rest,  and  a  spherlceuL  portion  of  the  fluid  is  suddenly  anni¬ 
hilated;  it  is  required  to  find  the  instantaneous  alteration  of  pressure 
at  any  point  of  the  mass,  and  the  time  in  which  the  cavity  will  be  filled 
iq>,  the  pressure  at  an  infinite  distance  supposed  to  remain  constant." 
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Raylel^  Integrates  the  equation  of  continuity  to  obtain  the  velocity 
distribution  for  the  flow*  He  derives  the  following  expression  for  the 
velocity,  ^2^ 

-  (1-1) 


u 


idxere  dots  Indicate  differentiation  with  respect  to  time,  R  Is  the  cavity 
radius,  and  r  Is  the  radial  coordinate*  For  the  motion  of  the  cavity 
wall,  with  the  use  of  the  principle  of  conservation  of  total  energy,  he 
obtains 


R« 


-Po 
P  , 


(1*2) 


where  p^  Is  the  pressure  In  the  fluid  at  rest  and  p  Is  the  density* 
Lord  Raylel^  also  derives  for  the  pressure  that 


^o  4' 


(1-5) 


Prom  Equation  1*2  R  =  U  =  -  oo  when  R  =  0*  Equation  1*3  shows 

that  when  R  «=  0,  the  pressure  p  Is  infinite  for  finite  values  of  r* 
Therefore,  the  energy  distribution  is  infinite  and  singular  at  that  in¬ 
stant*  Raylelfdi  also  shows  that  if  the  original  cavity  radius  Is  R^,  the 


collapse  time  is  proportional  to 


(l.U) 


and  is  therefore  always  finite*  Equation  1*4  follows  directly  from 
dimensional  considerations* 

Raylel^  realized  that  compressibility  should  be  considered  in  the 
later  stages  of  hlgdi  speed  collapse,  and  referred  to  work'by  Cook  who 
.  took  account  of  the  compressibility  of  the  cavity  contents  (he  assumed 
that  the  cavity  was  not  empty  and  that  the  contents  obeyed  Boyle's  law) 
but  not  of  the  surrounding  liquid*  Cook  examined  the  mode  of  collapse 
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when  the  cavity  strikes  a  rigid  sphere*  His  approximation  yields  plausible 
pressures  (68  tons  per  square  inch  in  a  pextlcular  example)  and  predicts 
that  the  collapse  velocity  will  be  zero  before  the  cavity  reaches  its 
mlnimimi  radius^  i.e*^  that  the  collapse  time  is  infinite.  However^  since 
his  work  does  not  account  for  the  compressibility  of  the  surrounding  medla^ 
it  is  not  satisfactory  from  a  theoretical  point  of  view. 

Investigations  concerning  the  roles  of  surface  tension  and  viscosity 

have  been  carried  out  by  Pcrit8ky(8)  and  also  by  Shu(9).  Poritsky  discusses 

the  collapse  or  growth  of  a  sj^erlcal  cavity  in  an  incompressible  viscous 

fluid  and  he  points  out  and  resolves  the  following  interesting  paradox  con~ 

cemlng  the  role  of  viscosity.  According  to  the  equation  of  continuity 

for  the  motion  of  an  incompressible  fluid  with  spherical  symmetry,  the 

velocity  ^  is  given  by  ^ 

V  .  (  )?■  ,  (1.5) 

r 


where  c(t)  is  an  admissible  time  function  and  is  a  unit  vector  along 

the  radius  r.  Hence  we  can  define  a  velocity  potential  cf>  as  follows. 


so  that 


^  =  V(f>  , 
<p  .  -clii 


(1-6) 

(1-7) 


Since 


0, 

V^/ir 


v^(v  (p)=  v(v^(p)^o 


(1-8) 


so  that  the  viscosity  term  in  the  Mavier  Stokes  equations 

grad  p  + 


Dv 


Dt 


=  F 


(1-9) 


vanishes  and  they  reduce  to  the  Eulerian  equations  of  motion 
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Poritsky  resolves  the  paradox  by  noting  that  while  it  is  true  that 
the  effect  of  viscosity  vanishes  in  the  equations  of  motion,  so  that  the 
resultant  viscous  stress  per  unit  volume  at  any  point  in  the  fluid  domain 
vanishes,  this  is  not  necessarily  the  case  with  the  stresses  themselves. 
At  any  point  the  three  principle  stresses  p^  and  strain  rates  ej^  are 
given  by  p 

Pi  =  “  P  -  3  »a(  +  eg  +  e^)  +  2tie±,  (l-lO) 

(i  =  1,2,3).  But  for  an  incompressible  fluid, 

e^  +  Cg  +  e^  =  0,  (1*11) 


and  hence  we  may  calculate  the  pressure  at  the  cavity  wall  from  Equation 
1.10.  We  may  not  calculate  the  pressure  from  the  pressure  in  the  fluid 
at  the  cavity  wall,  and  if  p  is  this  fluid  pressure,  then  for  the  pressure 
Po  at  the  cavity  we  have 

Po  «  P  -  2|ie.  (1.12) 


Here  p.  is  the  coefficient  of  viscosity  and  e  is  one  of  the  ccanponents 
ei-  As  Poritsky  emphasized,  viscosity  only  enters  the  problem  in  this 
boundary  condition. 

Poritsky  is  mainly  concerned  with  the  kinematics  of  the  collapse  and 
he  does  not  calculate  the  velocity  or  pressure  fields.  It  is  apparent  that 
if  compressibility  is  accounted  for,  then  the  role  of  viscosity  remains 
in  the  Navler-Stokes  equations.  Also,  we  no  longer  have  the  condition 
imposed  by  isochorlc  motion,  i.e..  Equation  1.11. 


Poritsky  finds  that 


Po  -  Pq> 


P 


R 

R  > 


(1.15) 


where  Po  is  given  by  Equation  1.12  and  Poo  is  the  pressure  in  the  fluid 
at  rest.  The  modifications  of  Equation  1.13  to  admit  the  effect  of  surface 
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tension  is 


R  R  + 


1  r2 

2  ® 


R 

T  I  • 


(1.14) 


where  CT  is  the  surface  tension  constant. 

He  finds  that  the  collapse  time  is  infinite  for  some  values  of  the 
coefficient  of  viscosity  and  Shu  establishes  rigorously  (with  modi- 
fications  of  the  Poincare-Bendixon  theory  of  non-linear  differential 
equations)  that  the  time  of  collapse  is  infinite  if  a  non-dimensional 
viscosity  c  is  greater  than  a  critical  value  and  is  finite 

otherwise.  Poritsky  shows  by  means  of  a  numerical  example  that  if  surface 
tension  is  accounted  for,  then  the  collapse  time  is  finite  and  Shu  claims 
that  this  can  be  shown  to  always  be  true.  Poritsky  also  shows  that  surface 
tension  speeds  up  the  collapse.  It  should  be  noted  that  Rayleigh's  theory 
always  predicts  finite  collapse  times. 

Zwick  and  Ple8set(l0)  investigated  the  dynamics  of  small  vapor  bubbles 
in  an  incompressible,  Inviscid  liquid.  Their  main  conclusions  may  be  sum¬ 
marized  as  follows.  The  vapor  pressure  at  the  cavity  wall  is  determined 
by  the  temperature  there.  Due  to  the  latent  heat  required  for  evaporation, 
a  change  in  bubble  size  will  lead  to  heat  transfer  across  the  bubble  wall 
causing  the  surrounding  liquid  to  heat  when  the  bubble  gets  smaller  and 
to  cool  during  bubble  growth.  Heating  the  liquid  causes  an  increase  in 
the  vapor  pressure  and  hence  slows  down  the  collapse. 

Tbe  results  obtained  by  Zwick  and  Plesset  show  that,  even  thou^  the 
temperature  of  the  bubble  wall  increases  rapidly  during  the  later  stages 
of  collapse,  the  motion  of  the  cavity  wall  agrees  very  closely  with  that 
predicted  by  Raylei^.  Hunter(ll)  concludes  from  these  results  that  the 
final  collapse  may  take  place  too  rapidly  for  the  vapor  to  condense  and 
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that  the  vapor  cavity  may  collapse  as  an  empty  cavity  until  the  vapor 
pressure  becomes  sufficiently  large  to  cactse  the  cavity  va]JL  to  rebound* 
Plesset(l2)  has  derived  the  equation  of  motion  for  the  bubble  radius 


R  in  an  inviscld,  incompressible  fluid*  He  finds  that 

E  R  *  I  S®  . 


(1*15) 


where  p  is  the  liquid  density,  Po  is  the  external  pressure  at  r  =  , 

and  p(R)  1b  tho  presanre-in  the  liquid  at  the  bubble  boundary  and  is 
given  by 

.  p(R)  =  Pv(T)  -  2cr/R*  (I*l6) 

Here  Pv(T)  is  the  equilibrium  vapor  pressure  corresponding  to  temperature  T 
at  the  bubble  boundary* 

The  effect  of  compressibility  on  the  collapse  of  a  spherical  bubble 
has  been  considered  by  Gllmore(l5),  Brand(l9)  and  Hunter(ll)*  Brand(l9) 
integrates  the  Eulerlan  equations  of  motion  numerically  using  the  method 
of  characteristics  and  compares  results  with  those  obtained  by  solving 
the  Lagranglan  equations  in  series*  A  Talt  equation  of  state  is  used, 
and  viscosity  and  surface  tension  are  neglected*  Hunter  also  carries  out 
a  calculation  by  the  method  of  characteristics  but  continues  up  to  the  in¬ 
stant  of  final  collapse,  while  Brand's  results  stop  short  of  this* 

Hunter  noted  from  his  numerical  results  that  the  motion  in  the  neigh¬ 
borhood  of  the  collapse  is  self  similar,  to  a  hi^  approximation.  However, 
in  formulating  his  similarity  theory.  Hunter  is  forced  to  approximate  vacuum 
ccmdltlons  at  the  cavity  wall* 

In  Hunter's  theory  there  is,  in  contrast  to  Raylei^'s  theory,  no 
singular  energy  distribution  at  the  time  of  collapse*  lEhe  similarity 
solution  predicts  an  infinite  velocity  at  the  time  of  collapse  but  gives 
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finite  pressures  for  all  radial  distances  at  this  time.  Hunter  shows  that 
compressibility  retards  the  collapse  process. 

The  collapse  of  a  cavity  in  an  incompressible  fluid  has  been  shown  to 
be  unstable  by  Birkhoff  (l4). 

The  present  paper  is  essentially  an  extension  of  the  work  of  Hunter. 

The  existence  of  a  self-similar  hypothesis  is  shown  to  be  implied  by 
dimensional  arguments.  Hunter’s  zero  order  boundary  condition,  that  the 
sound  speed  is  zero  at  the  cavity  wall,  is  corrected  by  use  of  a  per¬ 
turbation  scheme  which  is  developed  in  the  text.  Also,  several  of  Hunter's 
results  which  rely  on  energy  considerations  are  more  suitably  obtained  here. 
It  is  also  shown  that  the  cavity  wall  velocity  is  decreased  when  account 
is  taken  of  the  non-zero  sound  speed  at  the  cavity  wall.  Tiiis  collapse 
velocity  correction  is  obtained  without  numerical  integration  of  the 
equations  governing  the  first  approximation  to  the  non-self-similar  motion. 
Since  this  result  is  the  main  object  of  the  paper,  the  perturbation  equations 
are  not  integrated  here  although  the  procedure  for  doing  this  is  described 


in  full 
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2.  FOBMULATIOK  OF  THE  FROBLEM 

(1)  Similarity  and  Dlmensloiml  Techniques,  of  the 

Primitive  Eguatlona. 

Hie  symmetrical  collapse  of  a  cavity  la  water  Is  an  exampLeof  un¬ 
steady  motion  of  a  compressible  fluid  with  spherical  symmetry.  Under 
certain  conditions,  as  shown  by  Taylor(20),  Sedov(l5),  Stanyukovich(l6) 
Guderley(21)  £ind  other  authors,  the  problem  may  be  investigated  by 
similarity  techniques.  We  shall  discuss  the  concept  of  self-similar 
motion,  and  then  investigate  the  possible  self-similar  characteristics 
of  the  collapse  problem. 

A  motion  is  said  to  be  self-similar  if  the  spatial  distribution  of 
the  flow  parameters  at  a  certain  instant  is  identical  with  that  at  any 
other  time,  apart  from  a  change  in  scale. 

The  unsteady  motion  of  a  gas  with  plane,  cylindrical  or  spherical 
symmetry,  obeying  a  certain  type  of  equation  of  state,  can  be  shown  to  be 
self-similar  provided  that  the  number  of  constants  which  arise  with 
independent  dimensions  does  not  exceed  two.  For  completeness  we  shall 
include  the  discussion  of  this  point  given  by  Sedov.  Much  of  this  is 
relevant  to  Hunter's  analysis  of  the  cavity  collapse  problem. 

Let  us  analyse  the  dependent  variables  and  fundamental  parameters 
arising  in  unsteady  motion  in  one  space  coordinates.  In  the  Eulerian 
formulation,  the  physicsQ.  variables  may  be  teOcen  as  the  velocity  v, 
the  density  p,  and  the  pressure  p.  The  characteristic  parameters  are 
the  radial  distance  r,  the  time  t,  the  dimensional  constants  which 
enter  into  the  problem,  and  the  boundary  and  initial  conditions. 
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since  the  density  and  pressure  have  dimensions  irtiich  contain  the  mass, 
at  leeust  one  constant  a,  must  have  mass  occurring  in  its  dimensions.  We 
can  assume  that  the  dimensions  of  a  are 

[aj  =  ML^.  (2.1) 

nhen  we  may  write  for  the  velocity,  density,  and  pressure 


V 


^  “  k+1.  s+2 

r  z 


(2.2) 


^ere  V,  R  and  P  are  arbitrary  can  depend  on  non-dimensional  com¬ 
binations  of  r,  t  and  other  parameters. 

But  if  another  constant  b  entering  the  problem  has  dimensions 
independent  of  a,  and  all  other  constants  have  dimensions  depending  on 
a  and  b  only,  then  we  may  write 

[hj  ^  lV.  (2.5) 


Now  the  number  of  independent  variables  that  can  be  formed  by  combination 
of  a  and  b  is  reduced  to  one,  and  the  variables  v,  p,  p  will  depend  on 


this  one  non-dimensional  variable.  Moreover,  this  variable  will  be 


(2.4) 


and  if  m  #o,  this  similarity  variable  is  eq.uivalent  to  a  variable 

A  =  -1// -j-  ^  -  i  •  (2.5) 

b  t 

Hence  we  have  shown  that  if  there  are  only  two  dimensional  constants  with 
independent  dimensions,  then  the  motion  is  self-similar.  We  now  investigate 
the  eq^uations  governing  the  collapse  problem  and  detei.*mine  the  number  of 
constants  with  independent  dimensions.  * 


(2)  The  Role  of  the  Ener^ 


It  is  useful  to  clarify  the  role  of  the  energy  in  the  present  problem. 


Hie  total  energy  of  the  flaw  is  given  by 

QP 


U^J^TTr^clr^ 


(2.6) 


where  (3  is  the  internal  energy  and  u  is  the  velocity.  It  will  be  seen 
later  that  the  similarity  solution  to  the  present  problem  is  valid  within 
a  small  spherical  domain  (approximately  one  centimeter  in  diameter).  Hence 
energy  considerations  axe  strictly  inapplicable  in  determining  a  lower 
bound  for  the  similarity  parameter,  since  the  functions  u  and  given 
by  the  similarity  solution  cannot  be  used  to  compute  the  relation  given  by  Eq.. 
2.6  .  But  if  the  velocity  distribution  and  expressions  for  are  derived 
by  Integrating  the  governing  partial  differential  equations  exactly,  then 
we  can  impose  the  condition  that  energy  integral  (Equation  2.6)  should 
converge  uniformly. 


(5)  The  Ifodel 

-  • 

We  now  describe  the  model  of  collapse  and  the  dimensional  consequences 
of  this  model.  Following  Hunter  (see  footnote,  ibid.  pg.  246)  we  shall 
consider  a  cavity  initially  of  infinite  size  which  has  been  collapsing  for 
an  Infinite  amount  of  time.  As  he  points  out,  the  motiviation  for  this 
model  coaes  from  the  incompressible  treatment.  For  from  conservation  of 
total  energy  (using  Rayleigh’s  incompressible  treatment)  we  have 
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and  for  Ro»R^  we  see  that 

r2  «  2£^o -  ^  ^2 

5pR^ 

nierefore  the  motion  does  not  depend  on  the  scales  po  and  Ro  separately 

3 

hut  on  their  combination  and  since  is  eq.ual  to  3E  ,  •where  E  is 

tlT 

the  initial  energy  of  the  system,  and  if  we  let  Po  ->  o  and  Ro  -^oo  in 
such  a  way  that  E  remains  constant,  'then  the  flow  is  formally  determined 
by  "the  parameter  E  alone.  We  are  now  in  a  position  -to  formulate  the 
problem. 


(4)  Equations  of  Motion,  Boundary  Conditions,  Initial  Conditions 


Hhe  momentum  equation  and  "the  continuity  equation  corresixsndlng 
to  unsteady  motfon  of  a  gas  with  spherical  symmetry  are  respectively: 


-i.  ^.0 

(2.8) 

^  ■ 

(2.9) 

where  thft  primes  are  used  for  dimensional  quantities. 

The  Tait  equation  of  state  for  water  (22)  liay  bfi  Written  in  -the  form 

£4^  =  (  )*  .  (2.10) 


Here  B  and  Pq  are  slowly  varying  functions  of  entropy  and  they  may  be 
assumed  to  be  cons-tantj  if  is  -taken  equal  to  7  in  the  future  calculations. 
Qhe  velocity  of  sound  c'  is  given  by 


(2.11) 
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and  if  c’  «  Co  when  p*  =  Po>  then 

PpCo  *  ite* 


(2.12) 


Following  Hunter  we  shall  use  the  three  dimensional  scales  E>  po^  Cq 
to  express  all  variables  in  dimensionless  form.  We  then  have  (where 
variables  without  dashes  are  non-dimensional) 


d*  =  Mfio,  c*  =  c  Co>  P' 


E 


StTPo 


-  ^ 


(2.15) 


where  the  unprimed  quantities  are  dimensionless. 

73ien^  using  Equations  2.10,  2.11,  2.12,  and  2.1^,  Equations  2.8  and 
2.9  may  be  written: 


9u 


+  U 


du 

3? 


'St 


u 


dc^ 

ST 


1  dc^ 

y:3“ 

+  (i<-l)c2  ( 


0 


8u 


)  =  0 


(2.14) 

(2.15) 


The  boundary  conditions  are 

c  ■  1  and  h  =  R  at  the  cavity  wall,  r  =  R  (2.16) 

and  *  • 

c  ->1  as  r  00  (2.17) 

We  observe  that  in  the  incompressible  treatment  condition  (Equation  2*16) 
would  determine  the  motion.  The  effect  of  compressibility  is  felt  through 
the  initial  conditions  and  hence  the  spatial  variables  must  have  an 
assigned  distribution  at  a  given  time.  Hunter  uses  the  results  of 
Raylei^'s  incompressible  treatment  to  derive  initial  conditions.  From 
Equations  1.1,  1.2,  and  2.7  it  follows  that 


li^. 

5r 

These  initial  conditions  are  specified  as  t -oo,  vhen  the  fluid  can  • 

•*'k 

he  regarded  as  inco^ressihle*  Hence  the  realistic  conditions  in  the  final 

part  of  the  collapse  have  been  preserved  and  the  Equations  2.14,  2.15,  2.l6, 

2.17  and  2.18  completely  specify  the  flow. 

The  governing  Equations  2.8  and  2.9>  together  with  the  associated 

boundary  and  Initial  conditions  may  be  replaced  by  the  equivalent  system 

Equations  2.15>  2.l4,  2.15>  2.l6,  2.17  and  2.l8  .  We  then  see  that  three 

constants  with  Independent  dimensions  arise,  namely,  Cq,  po,  and  E.  To 

satisfy  the  requirements  of  similarity,  one  of  these  must  be  eliminated. 

Hunter  puts  equal  to  zero,  which  means  that  vacuum  conditions  apply  at 

the  cavity  wall.  Hunter's  Justification  for  this  step  is  that,  although 

the  dimensionless  sound  speed  at  the  cavity  wall  is  in  fact  equal  to  unity^ 

2 

in  the  region  where  similarity  solution  applies  c  1. 

In  the  present  analysis  we  take  account  of  the  fixed  sound  speed  at  the 
cavity  wall;  (note  that  Equation  2.12  shows  that  oq  =  0  implies  B  »  0  and 
hence  a  perfect  gas  law  for  water)  our  development  will  be  a  perturbation 
of  Hunter's  similarity  solution  such  that  his  solution  is  the  zeroth 
approximation  of  our  perturbation  scheme. 
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5.  CAVITY  COLT-APSK  WTTJB  A  NON-ZERO  SOUMD  SPEED  AT  rSE  CAVITy  WALL 
(l)  Perturbation  TOieory 

In  the  present  problem  three  dimensional  constants  enter  so  that  the 
similarity  hypothesis  Is  violated-  Several  investigators  have  dealt  vith 
similar  situations,  particularly  in  explosion  problems-  Sedov(15)  discusses 
the  point  explosion  vhile  taking  counter  pressure  into  account-  A  similar 
technique  vas  employed  by  Sakural(l7)  in  the  problem  of  a  blast  vave  prop¬ 
agation  vith  counter  pressure-  The  analysis  given  here  is  similar  to  that 
of  Sakural- 

We  define  ^  2  .2 

U  -  Rjjf(x,y),  c  =  Rgg(x>y) 

vhere  x  =  -^  and  y  =  s  ^  I 

U,  is  the  fluid  velocity,  ^  is  the  velocity  of  the  cavity  wall  as  given  in 

Hunter's  treatment  (i-e-  from  self-similar  equations),  and  c  is  the  sound 

speed-  Here  x  is  the  same  parameter  used  by  Hunter  in  his  self-similar 

solution  and  y  is  our  perturbation  variable-  All  quantities  are  dlmenslon- 

2 

less.  We  shall  expand  f  and  g  in  powers  of  y  as  follows; 

f  =  f^°^  +  y^f^^^  +  y^f^^^  +  .  .  .  .  (5-3) 

s  =  +  .  .  .  .  (3.4) 

where  f^^^  and  g^^^,  i  =  0,  l,.--^a  so  are  assumed  to  be  functions  of 
X  alone - 

We  are  assuming  that  Hunter's  solution  is  a  valid  zeroth  approximation 
and  hence  we  accept  his  determination  of  a  similarity  index  n  defined  by 

=1 - ^  a  constant. 


(5-1) 

(3-2) 


(5-5) 
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Hunter  finds  that  n  has  the  value  and  a  coinplete  discussion  of  his 

method  may  he  found  in  the  appendix,  (it  is  noteworthy  that  the  use  of 
Equation  5.5  simplifies  the  analysis  to  follow;  namely^  the  derivation  of 
the  non-self  similar  equations  to  determine  f^^^,  for  i  =  1,  2,  n>) 

By  substitution  of  Equations  5.1  and  5*5  in  Equations  2.14  and  2.15, 
by  use  of  the  expressions 

"Sr  R  ^ 

S  _  X  R  ^  d 
3t  "  "  R  ■  R  "Sy 

we  obtain,  after  some  reduction 

-X  ^  +  2/f  -i(y^  +  f^  -  ^  -  0, 

(5-8) 

and  • 

8g  +  f^  -  2fg^ 

^  OT  X 

+  g(i'-l)  (  ^  +  '■jc  '  1^)  =  °  (5.9) 

Substituting  Equations  5*5  and  5*4  into  Equations  5*8  and  5*9  respectively 

2 

and  comparing  coefficients  of  unity  and  y  we  obtain  the  zero  and  first 
order  approximations.  The  equations  may  be  arranged  as  follows: 


.  r(»)g(»)(l^)  (3-10) 


.  g(o)*  (M,  .  ,(o)^,Cc)  [  (, 


(3.11) 
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First  oi:der  approximation 

.(1) 


Here 


xf-  8^°^]  ^ 

xf-  g^°^]  g^^^PjCx)  +  f^^^^l^(x). 

1  X  *  X 

*  ^  8^  X  ■"  -  §Si)  -  4], 

Pj  -  2^x  -  +  xf^^^ki^-1)  +  f(°^(2  (^^-l)  -^(i^-1) 

.k/)  f(^f(. 

w  3C  3C  3C  3C  y 

Pj^  = 

+  1  g(o)  ^>)r2_  .2^  zJL\ 

^  8  [x  x^-1)  xj. 


(2)  Boundary  Conditions 

As  ve  have  seen,  the  exact  Boundary  conditions  are 
u  =  R  at  r  =  R 


and 


c  =  latr  =  R  , 


Let  us  first  consider  the  zero  order  approximation  .  Then  we  have 

.(o), 


and 


tt.-  Rgf(x,0)  =  §^''°^(x), 
d^  =  (x,0)  =  R|g^°\x). 


(3.12) 
\  -  •  •  •  . 

(3.15) 


(3.14) 


(5-15) 

(5.16) 

(3.17) 

(3.18) 
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TOie  functions  f^°^  and  are  solutions  of  Equations  5*10  emd  J.U.  They 

are  precisely  the  functions  that  occur  In  Hunter's  treatment.  We  see  that 
ve  are  unable  to  satisfy  the  exact  boundary  condition  (Equation  3.I6)  In  the 
zero  order  solution  (Equation  3.I8),  and  hence  to  this  degree  of  approximation 
we  have  g^°\l)  ■  0.  Obviously  condition  (Equation  3.15)  Is  satisfied  by 
setting  f^°\l)  *  1. 

We  now  substitute  the  complete  expansions  (Equations  3* 3  and  3*^) 

Into  Equation  3.1  and  apply  the  boundary  conditions  given  exactly  by 
Equations  3*15  and  3.16,  namely. 


R 


and 


^  [f^°''(l)+  +  .  .  .]  =  R, 

H  ■  i  •  .  j  -  1 

B  ^  B 


(3-19) 

(3.20) 


therefore.  If  we  set 

f^^^(l)  =  for  1  =  1,2,... n,  where  p^^^ls  to  be  (3*2l) 

determined,  and 

g^^\l)  =  1,  g^^\l)  -  0  for  1  =  2,  3,  ...n;  (3-22) 

then  we  can  satisfy  the  boundary  conditions  exactly. 

Therefore,  the  boundary  conditions  for  the  similarity  Equations 


3.10  and  3.11  are; 


and 

g^°^(l)  -  0. 


For  the  first  perturbation  of  the  self  similar  solution  we  have 
f(^)(l)  - 
g^^^(l)  -  1. 


(3.23) 

(3.24) 

(3-25) 

(3.26) 
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We  shall  now  discuss  the  zero  order  equations,  closely  following 

Hunter's  development*  Since  Hunter  does  not  tabulate  the  functions 
fio),  g^^^  It  Is  necessary  to  Integrate  Equations  In  the 

following  discussion  we  shall,  also  show  that  these  solutions  can  be  fixed 

without  appeal  to  the  energy  considerations  used  by  Hunter*  Further, 

Hunter's  discussions  of  the  singularities  of  Equations  3*10  and  3*11  will 

be  clarified* 


4*  THE  SELF  SIMILAR  BftUATIOMS 

In  the  appendix  a  complete  discussion  of  the  following  Items 
Is  given: 

(a)  ®ie  transformation  theory  employed  by  Hunter  in  his 
solution  of  the  s&Lf  similar  equations* 

(b)  A  discussion  of  the  singular  points  of  the  transformed 
equations  * 

(c)  The  application  of  techniques  given  by  Stanyukovich(l6) 
which  yield  the  same  conclusions  as  those  obtained  by 
Hunter  but  which  do  not  require  the  use  of  energy  arguments* 

(d)  A  discussion  and  results  of  the  numerical  integration  of  the 
transformed  equations*  (Hunter  gives  the  solution  curves 
but  he  does  not  tabulate  the  functions  f^°^  and  g^°^*  Since 
they  occur  in  the  coefficients  of  the  non-self  similar  equations 
it  was  necessary  to  repeat  the  Integration)* 
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In  Figure  1  ve  give  the  results  of  the  integration  of  the  self  similar 
equations;  Table  1  is  a  tabulation  of  the  function  f^°\  It  is 

Important  to  note  that  the  parabola  (f^°^  -  locates  the  singular 

points  of  the  f^°)  g^°^  and  f^^^,  g^^^  equations. 
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5-  AMALYSIS  OF  IHE  NON-SELF -SmiLAR  EQUATIOMS 


The  non-self -similar  motion  is  described  by  Equations  3*12  and  3-13* 
Gibese  equations  are  linear  differentisJL  equations  vlth  non-constant 
coefficients  having  singular  points  at  the  same  location  as  the  non-linear 
Equations  3*10  and  3"11;  namely,  at  the  two  points  where  the  quantity 
g^°^  -  (x  -  vanishes.  Obviously  this  vanishes  at  the  cavity  wall 

X  ••  1,  and  the  numerical  integration  of  the  non-linear  equations  shows 
that  the  expression  also  vanishes  at  x  1.31*  The  following  expansions 
are  immediately  obtained  from  the  non-linear  Equations  3*10  and  3-11* 

Near  the  cavity  wall  x  =  1. 


f^°^  =  1  -  2.0292(x-l)  +  .  . 
g(°^  -  4.8066  (x-1)  +  .  .  . 
and  near  the  singular  point  x  »  1.51> 

f^°)  =  .58  -  .57066  (x-1. 51)  +  • 
g^°^  =  .86  +  .409225  (x-1. 51)  + 


(5.1) 


(5.2) 


When  the  expressions  (Equation  5*1)  are  used  in  Equations  3"12  and  5-15^ 
the  following  expansions  are  found  for  f^^^  and  g^^^;  here  the  coefficient 
B  is  arbitrary  and  must  be  determined  from  uniqueness  considerations. 

We  have 


g^^^  »  1  -  11.354(x-1)  +  (x-1)^/^  -  2.979(x-l)^^/^+.  •  .) 


and 


,(1)  I’-4.807  +  4.369(x-1)  +  9.176(x-l)^l  f  „ 
^  “  \  -58.T6-$—  (x-ir 

b(i.166(x-1)^/^-  6.454(x-1)’^/^  +...)!. 

r28.326  -  ll6.47(x-l)  \  ^ 

1-28.1^5 Vl30.3fe(x-1)/  11.559(x-l) 

+  B  (  (x-l)"^/^-  2.979(x-l)^^/^  +  ••• 


(5.5) 


(5-4) 
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(l)  Analysis  at  the  Cavity  Wall 

We  now  consider  the  boundary  conditions  at  the  Cavity  weill.  The 
velocity  Is  given  by 


and  at  the  cavity  v8lL1 
R 


or  with 


-  i  2  •  •)> 


where  ^  Is  the  correction  to  the  self  similar  theory, 

+  « 


Prom  Equation  5*^  we  have  at  x  =  1, 

=  -.96567, 

and  hence  Equation  5*6  becomes 

^(1)  .96567 

and  therefore  the  new  velocity  of  the  cavity  wall  Is  given  by 

j,(i)  t  .96567 

B  ^ 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


(5.9) 


The  conqpirlson  between  R  end  ^  Is  shown  In  Figure  2. 


(2)  Analysis  at  the  Singular  Point  x  °  1.51 

With  the  use  of  expressions  (Equation  5*2),  Equations  5.12  and  5.15 
have  the  following  form  In  the  neighborhood  of  the  singular  point  x  =  I.5I: 
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,(1) 


(x-1.51)  +  Cg  (x-1.51)^  +  Cj  (x-1.51)5  ^  (5.10) 

*[A  +^3(x-1.51)^J  +  [pj^  +  P2(x-1.5l)+P3(x-1.5l)Jf^^} 

fc^(x-1.51)  +  C2(x-1.51)^+  C5(x-1.51)^J 


Here 


q2(x' 

-1.31)J  g^^^+  1 

jjl  +  r2(x-1.5l)+  rj(x. 

A- 

5.5178 

'll 

=  1.94511 

A  = 

9.072U 

=  -2.16913 

-24.39^ 

"•1 

=  .48o64 

Pi  “ 

.31514 

^2 

=  2.1183 

Pa  = 

6.8397 

^3 

=  -2.3324 

^5  = 

8.8287 

"1 

=  3.7933 

6.2369 

• 

"3 

=  2.4668 

At  x«1.51  expansions  for  and  given  by  Equations  5*10  and  5*11 


show  that 


g(l)^(x-l.5l)°^ 


(5.12) 

(5.15) 


where  and  Og  have  positive  non- integral  values.  In  fact  the  second 
order  equation  which  can  be  derived  from  Equations  5.10  and  5.11  has 
an  indlclal  equation  with  roots  2*79  and  .57^*  Therefore,  apart  from  the 
solution  =  0,  there  are  no  regular  integrals  at  the  singular 

point  X  *  1.51.  It  is  also  important  to  note  that  if  one  function  at  x  = 
1.51  is  specified  the  other  is  fixed,  since  for  example. 


2k. 


^1  +  PgCx-l.Jl)  +  P3(x-1.51)^^^^ 

-^3^(x-1.51)  +  c^ix-1.3lf  +  c^(x-1.3l)^} 


(5.14) 


and  there  exists  a  similar  expression  for  in  terms  of  and 


f(l) 

*  X. 


Hence  we  have  the  following  situation.  lEhere  etre  two  singular  points 
in  the  field.  Uie  expansions >  Equation  5.3  and  5*4,  describe  the  be¬ 
havior  of  and  g^^^  near  the  cavity  wall  singularity  and  eJ.so  f^^^  = 

-.96367  and  g^^^*  1.0  there,  but  there  is  an  arbitrary  constant  B  in 
the  expansions.  At  the  singular  point  x  =  1*51  there  are  an  infinity  of 
irregular  Integrals. 

Let  us  consider  the  apymptotic  behavior  of  solutions  of  the  per¬ 
turbation  equations.  For  large  values  of  x,  the  P's  in  Equations  3.12 
and  3.13  become 

p^.~  2Mx"*^ 

Pg - (2ffM+  1.3M)x"^*^ 

Pj^  2jfx+  (20-7M)x"^'^ 

-x~^'^(2jfM+  M) 


where  o  <  M  i  1  and  in  fact  M  o  as  x  Increases,  therefore,  from 
Equations  3*12  and  3*13  we  see  that  the  asymptotic  behavior  of  f^^^ 
and  g^^^  is  defined  by 


-  ^  MC^x~^ 

g<^Uc2X-^-602 


(5.16) 


where  c^^  and  Cg  are  arbitrary  constants. 
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Equation  5*l6  shovs  that  eJJ.  Integral  eurves  define  perturbation 
variables  which  tend  to  zero  at  infinity  so  that  the  similarity  solution 
iU  approached  as  the  distance  from  the  cavity  wall  is  increased* 

To  find  the  vaJLue  of  the  constant  B  and  the  pair  of  Integral  curves 
throu^  X  ■  1.51  a  further  condition  between  the  fluid  velocity  euid 
velocity  of  sound  must  be  satisfied  in  the  finite  part  of  the  field*  This 
can  be  derived  from  the  characteristics  solutions  of  Hunter  or  Brand. 

It  should  be  emihsisized  once  more  that  the  most  useful  product  of  the 
perturbation  scheme  is  the  correction  to  the  cavity  well  velocity,  ijuid 
this  is  found  without  Integrating  the  governing  equations  throu^butC  the 


field 
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APPENDIX 


Let  us  further  transform  according  to 

X  »  log(- J),  Y  =  -  fP,  Z  =  1‘^G.  (A-2) 

Then  we  may  arrange  the  transformed  equations  into  the  form 
dX:  dY:  dZ  =  (Y-l)^-Z  :  fY(Y-l)(nY-l)-  i  Z(9nY  +  n-lj| 

:2Z  C-nZ  +  7nY^  +  2Y(l-5n)  +  1 J  •  (A-5) 

V 

Let  us  consider  the  equation  for  ^  and  for  convenience  for  later  dis¬ 
cussion  define  y  =  Y-yo^  z  =  Z-Zq*  Then  from  (A-5)  we  have 


g , 

'’lo'  ^ 

where 

2  3 

P(y,z)  =  5nyz  +  y  (-n-1  +  5nyo)  +  ny  ; 

njifl 

Q(y^z)  =  yz(4-20n  +  28nyQ)  +  z^(-2n) 

+  y^(l4nZo)  +  zy^(l4nZo). 


(A-4) 

(A-5) 

(A-6^ 
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Also^ 


and 


1-n  ~ 

^lo  ‘  T  " 

aQl  “  -  2(n  +  l)yo  +  , 

^lo  “  2  +  yo(4-20n)  -  4nz„  +  l4nyo^, 
t>oi.  “  Zo('*-20n)  +  28nyoZo- 


Now  if  we  let 

y  =  vz, 

then  equation  (A-4)  becomes 


^  =  a*0  +  w  "  ^lo'^  ■  ^o.'^  ^YlCv/); 


dz 


where 


and 


2  3  3 

Vf  B  (-1-n  +  3nyQ)v  -  5iiv  +  nv^z  -  l^nzov"'^ 

^  2  3 

+  2nv  -(4-2C)n  +  28nyQ)v  -  l4nv''^z, 

-  l4nz  -  2n  +  (4-20n  +  28ny  )v  +  l4nv^z. 


But  for  z  »  0,  from  (A-9)  ve  see  that 

2 

va„.  +  a,  -  b,  V  -  b  ,v  =0. 

04  |o  lo  o.| 

If  we  let  be  a  root  of  (A-12)  and  further  define 
V  =  V  + 

then  equation  (A-9)  becomes 

z  *  Aq|  V  +  a  IQ  z  +  Ajjq  z  +  a  II  Vz  +  Aqp,.  + 

and  we  give  some  of  the  values  of  the  A's  as  follows; 

Ao,  ,  aQl  -  b,o  -  2V|  b„,^ 


^JO  +  ^0|  V, 


A,o 


(-1-n  +  3nyo)v,  -  Jnv,  -  l4nZoV,- 


2nv^  +  (4-20n  +  28ny^)V| 

^,o  +  ^o»  ^ 


(A-7) 

(A-8) 

(A-9) 

(A-10) 

(A-11) 

(A-12) 

(A-15) 

(A-14) 

(A-15) 


> 


(A-16) 
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^  __  nv^  -  l4nv^  -^4nz^Y,^-  2n  +(4-20n  +  28ny^)vjA|Q  ^ 

^  b/^  +  boi  V, 

p 

A„  -  2(-l-n  +3nyo)v,  -  3n  -  42nzQV, 

^10  +  ^OI  VI 

-  +2nv,  +  2(4-20n  +  28iiy^)v, 

^,o  +  ^ol  V. 

-  Ao,  (l4nzoV,^  -  2n  +(4-20n  +  28ny^)v, ) 

t|0  +  ^oi  V, 

-  A,o  L  (4-20n)zo+  28nyoZo] 


(A-17) 


(A-18) 


Equation  (A-14)  is  of  the  form  studied  by  Briot-Bouquet(l8)  and  ve  shall 
make  use  of  their  results  while  studying  the  solutions  of  Equation  (A-4) 
in  the  neighborhood  of  its  singular  points. 

(2)  Discussion  of  the  Singular  Points  of  the  ^  Equation 

•  - 

Prom  Equation  (A-5)  we  readily  find  the  aa&ilar  points  to  be  located 
in  the  (Y,Z)  plane  at 

(  (  Ion,  (A-19) 

2 

and  at  two  points  (E,D)  of  the  parabola  Z=(Y-1)  given  by  the  roots  of 

6nY^+  2(l-4n)Y  +  (l-n)  =  0  (A-20) 

Instead  of  using  energy  considerations  to  determine  which  of  the 
roots  of  Equation  (A-20)  is  appropriate  we  proceed  as  follows.  Since  the 
cavity  is  collapsing,  the  velocity  should  diminish  behind  the  front  of  the 
spherical  wave.  In  fact,  this  condition  is  necessary  for  the  existence 
of  the  similarity  solution  and  defines  the  domain  in  which  the  solution 


29. 


is  valid.  Analytically,  this  condition  states  that  in  a  nei^horhood 
behind  the  cavity,  ve  must  have 

dY>dX,  (A-21) 


and  this  should  hold  at  the  cavity  vail.  Due  to  the  incorrect  boundary 
condition  imposed  there,  the  inequality  (Equation  A-21)  is  undefined  at 
the  cavity  vail-  Hovever,  if  we  assume  Z  =  say  .05  there,  rather  than 
zero,  then  (Equation  A-21)  is  defined  and  shows  that  unless  we  choose 
the  larger  of  the  roots  of  (Equation  A-20)  we  will  have  n<.4  which  corres¬ 
ponds  to  the  incompressible  case. 

A  criticism  of  this  argument  is  that  it  leads  to  a  sound  speed  at  the 

cavity  weill  proportional  to  — However,  it  is  certainly  admissible 

^  ■ 

in  the  context  of  Hunter's  assertion  of  the  existence  of  a  self  similar 
motion  and  from  that  ix>int  of  view  it  is  believed  to  be  a  more  suitable 
argument  than  that  based  on  energy  considerations.  Moreover,  Hunter's 
numerical  Integration  shows  that  in  the  neighborhood  of  the  collapse,  R 
tends  to  a  constant  value. 

s 

We  recall  that  in  deriving  the  Brlot.Bcuquet  form  (Equation  A-l4) 
we  let  v,  be  a  root  of 

-b^jv^  +  v(a^  -  ^3,o)  +  a,o  =  0.  (A-22) 

Here  we  shall  show  conclusively  that  we  must,  ^at  the  point  E  defined 
by  Equation(A-20),  and  the  sonic  parabola  Z=(Y-1)^ 

have  the  root  of  Equation (a-12)  that  is  largest  in  absolute  value. 

First  we  consider  Equation  (a-20 ) .  From  this  we  see  that  the 
coordinate  Y  satisfies 


1 

¥ 


<  Y 


1 

n 


^  s  <  9 


(A-25) 


Hence 


(A-24) 
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cmd  If  we  transfer  the  location  of  the  singular  point  of  the  ^  equation 
to  the  origin  by  means  of  the  treuisformatlon 


(A-26) 

then,  at  the  origin,  the  Inequalities  (Equations  A-24  and  A-2^)  become 

(A-2T) 


1  <  yo<i 
¥  2 


1  <  Zq  <  9 

¥  °  T5 


Using  arguments  similar  to  those  establishing  the  dY  >  dX>  we  have^  at  E, 


dY>dZ.  (A-28) 

But  given  a  transformation  of  variables  of  the  form  (Equation  A-26),  It 
Is  easy  to  show  that  Implies  ^  >0.  For  If  a  function  t  =  p(q) 

undergoes  a  linear  transformation  which  maps 


and 


such  that 

and 


t  t 

q.  q 

(t  -  t) 

(q  -  q) 


are  of  the  same  sign,  then 


dt 


>  0 


d=E 


>  0 


and  conversely- 

Here  t  and  q  may  be  Identified  with  y  +  y^  and  z  +  respectively 
and  these  quantities  are  always  positive.  "Iherefore,  this  discussion 
shows  that  we  may  restrict  our  attention  to  Equation  (A-4)  which  defines 


Now  from  Equation  (A-12),  and  expressions  (A-7)  we  see  that 


2a 


tb 


*•01  i  ^  ^«0  "  ^lO  ^ol  ^ 


”  bio  > 

(A-29) 
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and  that 

V,  -  v,(n). 

Then^  since  y  «  vz. 


.  dv 


and  at  the  point  E  ve  have 

^  -  v,(n)>0. 

But  Equation  (A- 7)  and  the  fact  that  •4<n<l  shows  that  at  E, 


(A-3D) 


^  ®«>»  >  ^lo  >  ^oi  (A- 31) 

axe  eOJL  negative*  Hence  if  v,  is  to  be  positive,  then 

“o.  -  Ma„  b,,  -  «,.!>.,)<  0.  (A-32) 

But  with  the  range  assigned  to  the  variables  given  by  (A- 31)  it  is  easy 
to  show,  since  .4<n<l,  that 

(*o»  ^lo  “  ^o»)^*^*  (a- 33) 

Hence  unless  the  minus  sign  is  chosen,  v,  will  be  negative.  OMs  argument 
is  conclusive  whereas  Hunter  cannot  exclude  the  possibility  of  "exceptional 
ceises"  • 

There  are  three  pertinent  singular  points  of  the  W  equation,  and 

dZ 

when  we  transform  each  of  them  to  the  origin  we  have  the  following 
classification  of  the  singular  points  as  exhibited  by  the  nature  of  the 
characteristic  roots 
At  (l,0), 

2.667 
.  ..448 

and  hence  the  cavity  wall  is  a  saddle  point* 


52. 


At  (0,0) 

^1=1 

and  hence  the  origin  is  a  node. 

At  E(yo  =  .39^28,  Zq  =  .3669) 

>1  =  -.92755 
Xx«*  -.05855 

so  that  £  is  a  nodal  singularity  • 

Hunter's  paper  appears  to  have  a  misprint  concerning  the  drawing  of  the 
integral  curves  at  the  node  E  and  we  shall  include  the  correct  picture 
here. 


IHTEGRAL  CURVES  HEAR  NODE  E 
(5)  On  the  Integration  of  the  ^  Equation. 


In  order  to  integrate  the  equation 


dZ 


2Z(-nZ  +  7nY^  +  2Y(l-5n)  +  l) 
Y(Y-l)  (nY-1)  -  iz  (9nY  +  n-l). 


(A-3l^ 


it  is  necessary  to  investigate  the  behavior  in  the  nei^borhood  of  singular 
points.  We  get  appropriate  expansions  directly  from  the  theory  of  Briot- 
Bouquet.  They  show  that  if  in  Equation  (A-lU),  A^,  is  not  equal  to  a 
positive  integer,  then  the  Equation  (A-3^)  has  unique,  regular  expansions 


*  Note  that  we  Integrate  the  inverse  of  dY  since  Z  is  a  single  valued 

3z 


function  of  Y. 
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In  the  neighborhood  of  singular  points*  With  the  choice  of  the  root  v, 
at  the  node  £,  and  a  much  simpler  ansQ-ysis  at  the  other  singular  points^ 
it  Is  seen  that  is  never  a  positive  integer.  Then,  expanding  about 
a  singular  point  we  have 

Y  -  y©  =  (Z  -  2q)v,  +  c,(Z  -  Zq)^  +  C2(Z  -  Zq)^  +  .  .  (A-3S) 


where 


and 


Cl  = 


c 


2  = 


A|o 

l“Ao\  "  A,!  , 


Ao,c, 


+  A 


zo 


2  -  Aqi  -  A, I  . 


(A-3d) 


Then  Equation  (A-34)  i*e.,  the  W  equation,  has  the  following  expansions* 

dZ 

at  the  origin: 

Y  «=  .1Jj827Z  +  .008422Z^  +  .009882^  +  .  .  . 
at  the  point  E: 


Y  =  -.27701  +  2.166Z  -.969Z^  +  .  .  . 
at  the  point  (Y  =  1,  Z  =  0), 

Y  =  1  -  .hSlSSz  -.0964iz^  +  .  .  . 


(A-37) 


With  the  expansions  given  by  Equations  (A- 57)  can  integrate 
Equation  (A-  3^)-  iJiis  integration  was  carried  out  by  a  fourth  order  Runge- 
Kutta  technique.  The  results  of  this  integration  are  shown  in  Figure  V 
and  these  results  are  tabulated  in  Table  2. 
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TABLE  1 


TABULATION  OF  THE  FUMCTIONS  and 


X 

f(o) 

(0) 

X 

f(o) 

(0) 

g^ 

X 

f(o) 

1 

1 

0 

1.23 

.759 

.620 

1.46 

.615 

.835 

1.01 

.980 

.050 

l.g4 

.752 

.635 

1.47 

.610 

.840 

1.02 

.967 

.075 

1.25 

.745 

.650 

1.4a 

.665 

.845 

1.03 

.955 

.090 

1.26 

.737 

.665 

1.49 

.666 

.850 

i.o4 

•9^5 

.135 

1.27 

.729 

.680 

1.50 

.595 

.855 

1.05 

.935 

.200 

1.28 

.721 

.695 

1.51 

.590 

.86 

1.06 

•925 

.235 

1.29 

.715 

.710 

1.52 

.585 

.865 

1.07 

.910 

.265 

1.30 

.705 

.720 

1.53 

.580 

.870 

1.08 

.900 

.300 

1.31 

.700 

.730 

1.54 

.575 

.873 

1.09 

.892 

.342 

1.32 

.690 

.742 

1.55 

.575 

.875 

1.10 

.880 

.365 

1.35 

.685 

.750 

1.56 

.576 

.877 

1.11 

.875 

.385 

1.54 

.680 

.755 

1-57 

.565 

.879 

1.12 

.863 

.4l4o 

1.55 

.675 

.765 

1.58 

.560 

.881 

1.13 

.851 

.455 

1.56 

.670 

.775 

1.59 

.555 

.884 

1.14 

.840 

.470 

1.37 

.  665 

.785 

1.60 

.551 

.885 

1.15 

.830 

.490 

1.38 

.660 

.790 

1.61 

.548 

.886 

I.l6 

.820 

.510 

1.39 

.655 

.795 

1.62 

.544 

.887 

1.17 

.813 

.525 

l.4o 

.650 

.800 

1.63 

.542 

.888 

1.18 

.805 

.540 

1.4i 

.645 

.810 

1.64 

.538 

.889 

1.19 

.795 

.560 

1.42 

.640 

.815 

1.65 

.535 

.890 

1.20 

.785 

.575 

1.45 

.635 

.820 

1.66 

.531 

.891 

1.21 

.775 

.590 

i;44 

.625 

.825 

1.67 

.527 

.894 

1.22 

.767 

.607 

1.45 

.620 

.830 

1.68 

.523 

.895 
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X 

f(o) 

1.69 

•519 

•895 

1.70 

.515 

.894 

1.71 

•512 

.893 

1.72 

.510 

.892 

1.73 

.508 

.891 

1.74 

•505 

.890 

1.75 

.501 

.889 

1.76 

.498 

.888 

1-77 

•495 

.887 

1.78 

.491 

.886 

1.79 

.488 

.885 

1.80 

.485 

.884 

1.81 

.482 

.883 

1.82 

.480 

.882 

1.83 

.478 

.881 

1.84 

.475 

.880 

1.85 

.472 

.879 

1.86 

.468 

.878 

1.87 

.465 

.878 

1.88 

.460 

.877 

1.89 

.458 

• 

CD 

ON 

1.90 

.456 

.8750 

1-91 

.454 

.8725 

1.92 

.451 

.8700 

TABT.E  1  (  cont.  ) 


X 

,(o) 

(0) 
s'*  ’ 

1-93 

.449 

.8675 

1.94 

.446 

.8650 

1-95 

.443 

.8625 

1.96 

.440 

.8600 

1.97 

.438 

•8575 

1.98 

.435 

•8550 

1-99 

.432 

.8525 

2.00 

•  430 

.850 

2.01 

.428 

.848 

2.02 

.427 

.846 

2.03 

.425 

.844 

2.04 

K^ 

CVJ 

.842 

2.05 

.421 

.840 

2.06 

•  419 

.837 

2.07 

.417 

.854 

2.08 

.4i4 

.851 

2.09 

CVJ 

1— 1 

• 

.828 

2.10 

.410 

.825 

2.11 

.408 

.823 

2.12 

•  407 

.820 

2.13 

.405 

.818 

2.14 

CVJ 

0 

» 

.814 

2.15 

.400 

.810 

2.16 

.398 

.808 

2.17 

.397 

.805 

2.18 

•395 

.803 

2.19 

•  394 

.801 

2.20 

•393 

.798 

2.21 

•391 

.795 

2.22 

.388 

•792 

2.23 

•  385 

.789 

2.24 

•383 

•787 

2.25 

.380 

.785 

2.26 

•  378 

.783 

2:27 

.376 

.781 

2:28 

.374 

.778 

2.29 

.372 

.774 

2.30 

•  370 

.771 

2.31 

.368 

.768 

2.32 

.366 

.765 

2.33 

.364 

.762 

2.34 

.362 

•759 

2.35 

.360 

.756 

2.36 

•358 

•  752 

2.37 

.356 

•  749 

2.38 

.354 

•  746 

2.39 

.352 

•  743 

2.40 

•350 

•  740 

TABLE  1  (cont.) 
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X 

,(o) 

2.41 

.548 

.738 

2.42 

.546 

•736 

2.45 

.544 

.734 

2.44 

.345 

.732 

2.45 

.341 

.730 

2.46 

.540 

CD 

2.47 

.339 

.726 

2.4fi 

.337 

.724 

2.49 

.556 

CVI 

2.50 

.555 

.720 

5.00 

.135 

.310 

I 


20 


15 


10 
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FIG.  2  COMPARISON  BETWEEN 


UNCLASSIFIED 
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